Low density ferromagnetism in the Hubbard model 
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A single-band Hubbard model with nearest and next-nearest neighbour hopping is studied for d = 1, 
2, 3, using both analytical and numerical techniques. In one dimension, saturated ferromagnetism 
is found above a critical value of U for a band structure with two minima and for small and 
intermediate densities. This is an extension of a scenario recently proposed by Muller-Hartmann. 
For three dimensions and non-pathological band structures, it is proven that such a scenario does 
not work. 
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The Hubbard Hamiltonian (!]], which has recently at- 
tracted so much interest as a model for describing high- 
T c superconductors, was investigated already in the six- 
ties to face the problem of ferromagnetism in itinerant 
electron systems ^,|). Indeed, at a mean- field level, the 
Hubbard model seems to be a good starting point, as the 
Stoner criterion predicts a ferromagnetic ground state in 
a wide range of parameters. However, the inclusion of 
correlation effects makes the conditions for the appear- 
ance of ferromagnetism more stringent |2] . The existence 
of large-spin ground states in Hubbard type models has 
been rigorously established in three different but rather 
peculiar situations: Nagaoka ferromagnetism for infinite 
U and one hole in a half-filled band M, Lieb ferrimag- 
netism for half-filled bipartite lattices with sublattices 
containing a different number of sites and the flat 
(or nearly-flat) band ferromagnetism of Mielke || and 
Tasaki 

In a recent paper Miiller-Hartmann proposed an al- 
ternative route to ferromagnetism || . He considered the 
one-dimensional Hubbard model with both nearest- and 
next-nearest-neighbour hopping (ti-t^ Hubbard model). 
For infinite U and a band dispersion with two degener- 
ate minima he found saturated ferromagnetism for small 
electron densities. Very recently Penc and coworkers have 
studied a generalized version of Tasaki's model in lei and 
found ferromagnetism in a large region of parameters 
fiof . Their model reduces to that of Miiller-Hartmann 
in a certain limit to which, however, their considerations 
do not apply. 

In this paper we present numerical and variational re- 
sults for the Id t\-ti Hubbard model, which do confirm 
the existence of Miiller-Hartmann ferromagnetism. In- 
deed we find that this phase is not confined to infinite 
U and vanishingly small densities but occupies a sub- 
stantial portion of parameter space. It is then natural 



to ask whether this route to ferromagnetism works also 
for higher dimensions. We will show that in dimensions 
d > 2 the answer is negative. For low enough densities 
(and a non-pathological band structure), even in the pres- 
ence of several degenerate minima any state with finite 
magnetization is unstable. For the "marginal" dimen- 
sion d = 2 we cannot exclude ferromagnetism but rather 
restrict strongly the region of its stability. In fact, the 
existence of many degenerate minima is not enough to 
stabilize the fully polarized state: what really matters is 
a high density of states at the bottom of the band. 

We first consider the one-dimensional case. For the 
simple Hubbard model with only nearest-neighbour hop- 
ping, ferromagnetism has been discarded a long time ago 
by Lieb and Mattis, who showed that the ground state is 
a singlet pi. For the t\-t2 Hubbard chain 



H =- ^(iicf CT c i+iCT + t 2 cl a c i+ 2a + h.c.) + U^2 
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the Lieb-Mattis theorem is not applicable since when 
(2^0 the particles can pass each other and so cannot 
be ordered, which would be essential for the proof of the 
theorem. Actually, Mattis and Pena have shown in fj^| 
that, when t\ > and ti < 0, the model shows Nagaoka 
ferromagnetism for one hole in a half-filled band. Here 
we are interested in small and intermediate densities. 

We have performed numerical diagonalizations using 
the Lanczos method for chains up to L = 18 sites, for 
even electron numbers N < L (n = N/L is the electron 
density). The total spin of the ground state is determined 
both by measuring the spin-spin correlation function at 
q = and by calculating the lowest energy states for 
each S e subspace. The energy of the fully polarized state 
is independent of U and calculated analytically. When 
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t2 > —t%/4 (we always suppose t\ > 0) the band has 
only one minimum at k = and we do not find ferromag- 
netism at low density: for every value of U the ground 
state is a singlet. From now on we consider only the case 
t% < — ii/4, when the band developes two degenerate 
minima at momenta ±ko, with fco = arccos(ii/4i2). It is 
in this region that we should find, according to Miiller- 
Hartmann, ferromagnetism at low density and infinite U. 
We note that the lowest density that can be studied on a 
finite chain is 2/L. For this case of two electrons, one can 
easily show that the ground state is a triplet. However, 
this result can be also regarded as an example of Mielke's 
flat band ferromagnetism. Interestingly, for large enough 
U, we do find ferromagnetism also for more than two 
electrons, namely for a range of densities up to a criti- 
cal value n c . This is no longer flat band ferromagnetism 
but rather confirms the validity of Muller-Hartmann's 
construction beyond the limit n — > and U = oo. 

In our numerical diagonalizations we have never found 
partially polarized ground states: the phase transition is 
always between a singlet and the fully polarized state. 
Fig.l shows the behaviour of function of tz/ti. 

Data from chains of different lengths are plotted in the 
following way: if for instance with N = 4 electrons the 
ground state is ferromagnetic, while with N = 6 it is 
not, we draw an error bar between 4/L and 6/L. The 
solid curve represents the density at which the number 
of Fermi points for the fully polarized state changes from 
four to two. The curve seems to fit the data rather well 
and shortly, we are going to present an analytical argu- 
ment that this is the likely criterion for the phase bound- 
ary]^]. In any case, one can say that for densities above 
this curve Muller-Hartmann's construction is no longer 
applicable since it relies on the existence of two discon- 
nected pockets. In summary, our results indicate that 
this type of ferromagnetism is not restricted to very low 
densities, but it survives as long as the (fully polarized) 
Fermi sea has four Fermi points. 

We turn now to the question which value of U is re- 
quired to sustain ferromagnetism at a given density n. 
Since our numerical results did not show evidence for 
a partially polarized ground state, we may examine the 
stability of the fully polarized state \F) with respect to a 
single spin flip. We use the standard single-spin-flip trial 
state El 



(2) 



An up-spin electron is removed from one of the Fermi 
points hp (with energy ep) and a down-spin electron is 



*At least for not too large |ta|. In fact, for |ta|/ti 
system decouples into two simple Hubbard chains. 
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put to the bottom of one of the degenerate minima eo at 
wave vector fco. This yields a single-particle energy gain 
of A = e F — €q. The change of the total energy is 
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where Eq = E\ + E2 is the energy of the fully polar- 
ized state, while E\ and E2 are the contributions from 
the nearest-neighbour and next-nearest-neighbour hop- 
ping terms, respectively. AE has to be minimized with 
respect to the variational parameter rj. 

The polarized state becomes unstable when the kinetic 
energy gain outweighs the increased interaction energy. 
Fig. 2. shows, for the special value t<z = —t\ the vari- 
ational phase boundary U vax {n) (continuous curve) and 
the critical U c determined by the numerical diagonaliza- 
tion (dots). We found similar results for other values of 
t% (< — ti/4). We notice that Z7 V ar is a lower bound for 
the exact Uc. For small density n, U c decreases with de- 
creasing n, apparently tending to zero as n — » 0. The 
cusp of J7 var (71) occurs at precisely the critical density we 
have defined above and indicates that something drastic 
happens there. While for n > n c we see a sharply rising 
t/ var (n), our numerical data indicate that ferromagnetism 
is completely suppressed in this region. Since the varia- 
tional result provides a quickly increasing lower bound for 
U c , we take it as an indication that, in fact, U c = 00 for 
n > n c . In this sense, the variational approach supports 
the phase diagram shown in Fig. 1. 

We investigate now the stability of low-density ferro- 
magnetism for d > 2. We consider the Hubbard Hamilto- 
nian with arbitrary non-diagonal hopping terms tij and 
vanishing diagonal terms, ta = 0. However, we confine 
ourselves to non-pathological band structures, meaning 
quadratic dispersion about the minima. Since the most 
favourable situation for a fully polarized state is U = 00, 
we use the variational ansatz (0) for t] = 1. A straight- 
forward calculation gives 
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where A = ep 



to ~ n 2 / d . (In Id, the leading cor- 
rection would be 0(n 2 ).) We note that eo, the lowest 
eigenvalue of the hopping matrix, is necessarily negative, 
as the sum of the eigenvalues is equal to J2i Ui — 0. The 
first term in Eq. (||) is therefore positive, favouring ferro- 
magnetism, while the second term destabilizes the fully 
polarized state. 

We examine now the viability of the low-density route 
for dimensions d > 1 . For d > 3 the second term in (^) is 
the leading term at low density and the polarized state is 
always unstable for sufficiently low densities. (For d = 1 
the first term is dominating and the present argument 
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does not predict a destabilization of the fully polarized 
state.) For d = 2 the two terms are of the same order 
and must be compared more carefully. In this case, to 
leading order in n, A = n/p(eo), where p(eo) is the den- 
sity of states at the bottom of the band. Inserting this 
expression into Eq. (Q) and using the stability criterion 
AE > 0, we obtain the following necessary condition for 
the stability of the fully polarized state, 

2p(e )|e | > 1. (5) 

To be specific, let us consider the t\—t^ Hubbard model 
on a square lattice. When ti < —t%/2 the band has de- 
generate minima at the Brillouin zone boundary. It fol- 
lows from Eq. (|5|) that the polarized state cannot be the 
ground state outside the region — 0.20ii > t<i > — 0.65ti. 
Hence for t% < — 0.65ii the polarized state is unstable 
even though there are two minima]]]. 

Having excluded low-density saturated ferromagnetism 
for d > 2, we may wonder whether we can also say some- 
thing about partially polarized states. Unfortunately, en- 
ergy eigenstates with S < S ma x are not exactly known, 
so that we have to find a different way for judging the 
stability of ferromagnetic order. An essential ingredi- 
ent for the argument is the monotonous increase of the 
lowest energy eigenvalue as a function of U, within any 
given subspace. This holds for the lowest singlet energy 
E S (U) as well as for the lowest energy E pp (U) within a 
subspace with a given partial polarization, as depicted in 
Fig. 3. An additional important fact is that the Hartree- 
Fock state gives an upper bound for the the lowest sin- 
glet state Ehf{U) > E S (U). Therefore the energies for 
the partially polarized state and the lowest singlet state 
cannot cross at a smaller U than the value defined by 
E HF (U) = Epp(0) (see Fig. 3). But this value of U be- 
comes arbitrarily large as the density goes to zero. In 
fact, for d > 2 it diverges as n 2 / d ~ x . 

We should comment on the relationship between our 
description of the low-density route to ferromagnetism 
and other scenarios of ferromagnetism in Hubbard mod- 
els. The existence of the single-hole Nagaoka state is 
well-known jlj] but it seems clear that the fully polar- 
ized state of the Id t\—t% model at low to intermediate 
densities is in its physical nature different, and in the 
phase diagram disconnected, from the Nagaoka state. In 
the special case n = 2/L, our ferromagnetism is indistin- 
guishable from Mielke's [0] but we find that the ordered 



^The t-matrix expansion has been also used to study the 
possibility of low-density ferromagnetism in the 3d and 2d 
Hubbard models Jl^JiB) . Let us remark however that this 
expansion, even if physically justifiable, is nevertheless un- 
controlled. Moreover, the role of several equivalent minima 
was not discussed. For the case of a single minimum, our 
results are in agreement with those of H,f5| • 



state extends up to a ^/independent ) critical n c , and 
for finite concentrations n we certainly do not have flat- 
band ferromagnetism. 

We have already mentioned that the generalization of 
Tasaki's model |1 by Penc et al. pi includes Muller- 
Hartmann's model || as a special case. However, the 
understanding of ferromagnetism in Jl0[ is based on the 
perturbative treatment of ring exchange processes, and 
this scheme breaks down for the model considered by us. 
Our finding that a decisive role is played by the exis- 
tence of four Fermi points, seems to have no counterpart 
in their reasoning. The possible relationship between the 
ferromagnetic phase discussed in [ fiof and the one de- 
scribed by us, remains to be clarified. 

It is worthwhile to point out that the existence of 
several degenerate minima does not play an essential 
role in our reasoning for d > 2. The phenomenon of 
valley-degeneracy-assisted ferromagnetism is apparently 
restricted to Id systems. 

To conclude, we have found an extended region of fully 
polarized ferromagnetic states in the U-n plane for the 
Id t\— 12 model. The available evidence indicates that 
ferromagnetism exists whenever the Fermi sea consists of 
two disconnected parts. The ferromagnetism described 
by us is not a flat-band phenomenon, and it does not 
seem to be connected with the Nagaoka state, either. 
It is most unfortunate that the mechanism giving rise 
to the spin-polarized state in Id does not work for non- 
pathological situations in 3c? and has, at best, a rather 
restricted range of applicability in 2d. 
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Figure Captions 

Fig. 1. Phase diagram of the ti-t 2 Hubbard chain in 
the ti~n plane, as obtained by exact diagonalization for 
finite lengths L. Circles: L — 10, diamonds: L = 12, 
squares: L = 14, triangles: L = 16. The full line shows 
the density where the fully polarized Fermi sea splits into 
two pockets. 

Fig. 2. Phase diagram of the Hubbard chain 

in the n—U plane, for ty. = —t\. Circles show numerical 
results, while the full line gives a variational lower bound. 
The location of the cusp (n c = 2/3) corresponds to a 
point on the full line in Fig. 1. 

Fig. 3. Schematic view of the {/-dependence of the low- 
est energies for subspaces with different values of the total 
spin (s: singlet, pp: partially polarized) . The dashed line 
shows the non-magnetic Hartree-Fock energy (HF). 
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